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Abstract 

We summarize recent results on i) the use 
of u-pdf 's in the Monte-Carlo simulation of 
hadronic final states at high-energy colliders, 
and ii) attempts to characterize u-pdf 's with 
precision in terms of nonlocal operator ma- 
trix elements. 

1 Introduction 

At forthcoming high-energy colliders a large 
amount of events are characterized by com- 
plex final states with multiple hard scales, 
possibly far apart from each other. QCD 
methods to address multi-scale hadronic pro- 
cesses involve the use of parton distributions 
unintegrated in both longitudinal and trans- 
verse components of parton's momentum (u- 
pdf's). Classic examples are provided by Su- 
dakov physics [BE] and small-x physics [3l[4]. 
U-pdf's also govern exclusive processes [5] 
and spin [6]. They are required, in general, 
to explore detailed features of hadronic final 
states [7]. For these reasons, theoretical and 
phenomenological studies of u-pdf's are be- 
ing actively pursued. 

In the case of small x, u-pdf's can be in- 
troduced in a gauge-invariant manner using 
high-energy factorization [8j . This result was 
used early on both for Montc-Carlo simu- 
lations [9] of x ^ parton showers and 
for numerical rcsumniation programs |10j for 
In X corrections to QCD evolution equations. 
For structure function's evolution, methods 
are being developed ([ll]-[T5] and references 
therein) to match the k^-dependent, small- 
X dynamics with perturbative coUinear dy- 
namics. For the full simulation of exclu- 



sive components of hadronic final states, on 
the other hand, such matching is more com- 
plex (see discussions in [7] and [16]) and is 
not yet available. This will be critical for 
turning present event generators based on u- 
pdf's [17]-[23] into general-purpose Monte- 
Carlo tools. Importantly, for x — *■ u- 
pdf's may provide a useful framework also 
for discussing the approach to the saturation 
regime: see e.g. [531 123 for recent studies. 

In the general case, to characterize u-pdf's 
gauge-invariantly over the whole phase space 
is more difficult. Open questions on this is- 
sue are the subject of much current activity, 
see e.g. [26]- [40], with wide-ranging applica- 
tions from semi-inclusive processes to parton- 
shower algorithms to spin physics. 

In this report we concentrate on two top- 
ics. In Sec. 2 we consider x <C 1 and discuss 
recent results from k_L Monte-Carlos on the 
jet structure of small-x final states. In Sec. 3 
we describe ongoing progress towards general 
operator definitions for u-pdf's, focusing on 
results for the treatment of x — > 1 endpoint 
divergences. 

2 Jet final states from uninte- 
grated parton distributions 

There have recently been new calculations of 
hadronic jet final states from Monte-Carlo 
generators based on unintegrated parton dis- 
tributions. 

2.1 U-pdf's and shower Monte-Carlo 
generators 

The idea common to the Monte-Carlo event 
generators based on u-pdf's is to use factor- 
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ization at fixed k_L [5] in order to a) gener- 
ate tlie hard scattering event, including de- 
pendence on the initial transverse momen- 
tum, and b) couple this to the evolution of 
the initial state to simulate the gluon cas- 
cade. Event generators of this kind include 
SMALLX 0], CASCADE [ID], LDCMC [12], 
and the newer tools and [18]. Differ- 
ent generators employ different models for 
the evolution of the initial state, such as 
BFKL [4ll[42], CCFM [43l[44], LDC [45] evo- 
lution. With suitable constraints on the an- 
gular ordering of gluon emissions, these can 
all be set up to give the correct leading Inx 
behavior. A number of subleading contribu- 
tions is however also of importance for real- 
istic simulations of final states. The general 
approach of these Monte-Carlos and aspects 
of the differences among them are reviewed 
in [16]. Sec also [7] for further references. 

Although none of the above generators 
are nearly as developed as standard Monte- 
Carlos like PYTHIA [M] or HERWIG [17], 
they have the potential advantage of a better 
treatment of high-energy logarithmic effects. 
Implementing these effects in the shower can 
be relevant for the simulation of complex fi- 
nal states at LHC energies, see e.g. multi- 
jet studies in [48], and additional references 
in [7]. A further potential advantage is 
that Monte-Carlos with u-pdf's likely pro- 
vide a more natural framework [49j to sim- 
ulate the k^ distribution of the soft under- 
lying event |50| (minijets, soft hadrons) and, 
possibly, multiple interactions at very high 
energies. 

On the other hand, present Monte-Carlos 
of this kind do not automatically include con- 
tributions from coUinear radiation associated 
to X 1. They need to be corrected for this. 
This is done partially in present implemen- 
tations, but not yet in a systematic fashion. 
A recent example of a study in this direction 
is in [17] , based on the procedure of [ST] [52] . 

A related issue concerns the inclusion of 
quark contributions in the initial state. In 
present implementations these are either ne- 



glected or approximated by lowest-order per- 
turbative evolution [H [THl [H]. The k_L- 
dependent kernel that governs small-x sea- 
quark evolution to all orders is given in |53j . 
Universality properties of this kernel are em- 
phasized in [l2l HI] . The quark kernel is not 
yet implemented in current Monte-Carlos. 

Keeping in mind these limitations, Monte- 
Carlo event generators with u-pdf's can be 
used at present within the low-x domain to 
probe the physical picture of space-like x ^ 
parton showers. A selection of recent results 
from these event generators is given in the 
next subsections, going from more inclusive 
to less inclusive measurements. 

2.2 Inclusive cross sections 

Available k^ Monte-Carlos have been used 
to analyze data for DIS structure functions 
and for inclusive lepto- and hadro-production 
of jets, heavy flavors and prompt photons, 
see [7] [IS] . These inclusive analyses serve to 
test the overall consistency of the physical 
picture and to perform first determinations 
of the unintegrated gluon distribution. 

A consistent semi-quantitative picture is 
achieved with sensible results for the evolved 
unintegrated gluon distribution [THl [IHl Ell 
[SH [51 [H [55]- However, the gluon at low 
scales and low x is only poorly constrained 
by this approach [58] [59] . This is not unex- 
pected, as the summation of higher orders 
implied by the Monte-Carlo should reduce 
the sensitivity of predictions at high scales 
to the form of the input [60] . 

Examples of results for jet production in 
ep and pp collisions are reported in Figs. [1] 
and [2] Fig. [T] shows the description of the 
jet Et distribution data from HI [61] ob- 
tained with CASCADE using the recent u-pdf 
fits [52]. Fig. [Ushows the comparison of CDF 
jet data [63] with the Monte-Carlo [17]. 

Current uncertainties on the unintegrated 
gluon are large especially in the lowest x 
region [16j . Implications of this are stud- 
ied for heavy-quark structure functions at 
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Figure 1: Jet Et distribution at HERA from 
the k_L Monte-Carlo CASCADE and u-pdf 
fits[62], compared with HI data|61|. 
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Figure 2: Tcvatronjct spectra from|17| com- 
pared with CDF data[53]. 



terms, although highly subleading from the 
viewpoint of the soft hierarchy, appear to af- 
fect the spectrum at a level comparable to the 
current theoretical uncertainties estimated as 
in [67] . Further investigations of these ef- 
fects [221 EH ESI EQl [71] are warranted. An 
additional, special application is to central 
exclusive production, recently considered for 
Higgs [75] [73] and scalar x mesons [71] in the 
context of u-pdf's. 
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Figure 3: (top) Azimuth and (bottom) 
Bjorken-x dependences of di-jet distributions 
measured by Zeus [78] . 



HERA in [58] and for jet production at RHIC 
in [ST] . Note that first fits of both the x- and 
k_L-dependence of the unintegrated gluon are 
performed in [62j . 

A potentially important application of 
these results is to final states containing 
Higgs bosons at the LHC. The Higgs matrix 
element coupled to the unintegrated gluon 
distribution is computed in [M]. First re- 
sults from Monte-Carlo implementation are 
given in . Studies of the Higgs transverse- 
momentum spectrum from NLO + soft-gluon 
resummation |651 I66| indicate that small-x 



We conclude this subsection by observing 
that recently the small-x definition [5] of u- 
pdf's has also been used [75l|76] to treat non- 
leading corrections to jet production at large 
rapidity separations in hadron-hadron colli- 
sions [77], and Monte-Carlo event generators 
based on the evolution [JTJ [32] are also be- 
coming available for these processes. 

2.3 Multi-jet correlations 

More details of the parton k_L dynamics can 
be probed by examining jet correlations in 
final states containing multiple jets. 
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Recently the Zeus collaboration has pre- 
sented measurements of two-jet and three-jet 
distributions in DIS associated with low x, 
10^4 < X < 10~2 [78]. Zeus' jet definition 
is such that nonglobal logarithms from jet 
clustering are avoided [79] , as are double log- 
arithms from symmetric cuts [SO]- Fig. [51 [75] 
shows two-jet data versus the azimuthal sep- 
aration A(/) and versus Bjorken x, compared 
with next-to-leading-order perturbative re- 
sults [ST]. Notice the large variation from 
the order-a^ to the order-a^ result with de- 
creasing X and decreasing A(/). 




Figure 4: Angular jet correlations [5^ from 
CASCADE and HERWIG compared with Zeus 
data[78]: (top) di-jet cross section; (bottom) 
three-jet cross section. 

Ref. [60] analyzes the multi-jet final states 
in the framework of u-pdf 's. In particular, it 
studies the contribution to potentially large 
higher-order corrections arising from a size- 
able k_L in the initial state when several well- 
separated hard jets are produced. In this 
kinematics, effects from the x ^ 1 parton 
shower may be enhanced. These contribu- 
tions are important at small A0, when the 
jets arc not close to back-to-back configura- 



tions [52] . Fig. [3] shows results for the distri- 
bution in the azimuthal separation between 
the leading jets |60| . compared with the mea- 
surement [75] . Fig. [S] [SD] shows the angular 
distribution for the third jet, in the cases of 
small separations and large separations be- 
tween the leading jets. The kj^ Monte-Carlo 
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Figure 5: Cross section in the azimuthal an- 
gle between the hardest and the third jet [60], 
for (top) small and (bottom) large azimuthal 
separations between the leading jets. 

CASCADE gives a good description of the 
measurement, and gives large differences to 
the collinear-based parton shower as imple- 
mented in HERWIG, particularly in the re- 
gion where the azimuthal separations be- 
tween the leading jets are small. 

The X dependence of the gluon distribution 
used in the above results for the jet cross sec- 
tions is shown in Fig. [6] for various scales /i 
and fixed k^. 

This analysis indicates that theoretical un- 
certainties on jet correlations are sizeable at 
NLO [78] and that the kj^ Monte-Carlo gives 
an improved description of multi-jets [60j 
compared to standard parton showers. This 
analysis is for DIS. However, despite the 
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Figure 6: The x dependence of the unintc- 
grated gluon distribution. 

much lower energy at HERA, note that ow- 
ing to the large phase space available for jet 
production the HERA data may be just as 
relevant as the Tevatron for understanding 
the extrapolation to the LHC of initial-state 
radiation effects [7]. 

2.4 Forward region 

Transverse energy flow and production of jets 
and particles at forward rapidities are im- 
portant probes of the initial-state kj^ [55] , 
Monte-Carlo results for such observables, 
however, strongly depend on the details of 
the model used for evolution of the u-pdf's. 
This model-dependence is much more pro- 
nounced for forward-region observables than 
for the jet correlations discussed in the previ- 
ous subsection. A systematic understanding 
of u-pdf's evolution [M] will be relevant for 
more precise interpretations of forward mea- 
surements. See [40] for recent work on evolu- 
tion equations, including discussion of target 
fragmentation. We refer the reader to [3 [H] , 
and references therein, for discussion of cur- 
rent forward-region results. 

3 Towards precise characteriza- 
tions of u-pdf's 

The possibility to turn the event genera- 
tors discussed in Sec. 2 into general tools 
to describe hadronic final states over the 



whole phase space depends on theoretical 
progress about unintegrated parton distri- 
butions. This will involve precise charac- 
terizations of the distributions, factorization 
and evolution equations. In what follows we 
discuss recent progress in the understanding 
of operator matrix elements for u-pdf's and 
their lightcone limits. 

3.1 Gauge invariant matrix elements 

The relevance of consistent operator defini- 
tions for parton k^ distributions was em- 
phasized long ago, see e.g. O [85]. Ensur- 
ing gauge invariance at k^ ^ 0, however, is 
nontrivial. To this end the approach com- 
monly used is to generalize the coordinate- 
space matrix elements that serve to identify 
ordinary pdf's to the case of field operators 
at non-lightconc distances. E.g.. for quarks 
one has [ilMlIlIj 

fiy) = {P\^iy)V^in)j+Vo{n)iP{0)\P) . 

(1) 

Here tp are the quark fields evaluated at 
distance y = {0,y~,y±), y± is in general 
nonzero, and V are eikonal-line operators 
in direction n required to make the matrix 
clement gauge- invariant. The unintegrated 
quark distribution is obtained from the dou- 
ble Fourier transform in y~ and y±^ of /. An 
extra gauge link at infinity |86| is to be taken 
into account in the case of physical gauge. 

There are subtleties, however, to using 
Eq. ([T]) beyond tree level. As realized early 
on |85| . parton distributions at fixed kj_ 
are no longer protected by the KLN mech- 
anism [55] against uncancelled lightcone di- 
vergences near the x = 1 endpoint ^87j i89j . It 
is only after supplying the above matrix ele- 
ment with a regularization prescription that 
the distribution is well-defined. 

Importantly, similar to what observed 
in [90] for the case of the Sudakov form fac- 
tor, the choice of a particular regularization 
method for the lightcone divergences also af- 
fects the distributions integrated over kj^ and 
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the ultraviolet subtractions. A further set of 
questions being studied [571 EH IHIl are in 
fact associated with the coefficient functions 
that govern the expansion of unintegrated 
distributions in terms of ordinary ones, and 
the relation of the integral of u-pdf's with 
ordinary distributions. 

The above issues can be analyzed by ex- 
plicit calculation at one loop [27l [93] . Ex- 
pansion in powers of of the result for the 
coordinate-space matrix element at this or- 
der (Fig. [7]) yields 



We consider these methods briefly in the next 
subsection. 
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where /i is the dimensional-regularization 
scale and p is an infrared mass regulator. The 
lightcone singularity v —^ I, corresponding to 
the exclusive boundary x = 1, cancels for or- 
dinary pdf 's (first term in the right hand side 
of Eq. ([2])) but it is present, even at d 4 
and finite p, in subsequent terms. This im- 
plies that, using the matrix element ([T|), the 
1/(1 — 2;) factors from real emission probabil- 
ities do not combine with virtual corrections 
to give 1/(1 — .t)_|_ distributions, but leave 
uncancelled divergences at fixed k_L. 



(a) (b) 
Figure 7: One-loop contributions to lightcone 
divergences in the matrix element ([T|). 

Different methods have been employed in 
the literature to provide u-pdf's with consis- 
tent regularization for the cndpoint region. 



3.2 Cut-off vs. subtractive method 

A possible approach to treat the endpoint is 
to implement a cut-off by taking the eikonal 
line n to be non-lightlike, as in the early 
works [85] and [94] . The cut-off in x at fixed 
k^ is of order 1 — x > kj_/ where 77 = 
{p-nY' /"n? and p is the incoming momentum. 
Evolution equations in the cut-off parame- 
ter 77 are investigated in [3S1 [5S1 [^3 [HS] . A 
thorough analysis of factorization at leading 
order has recently been given in [26j based 
on this approach. 

A potential drawback of this approach is 
that cut-off regularization is not very well- 
suited for applications beyond the leading or- 
,„sder. Also, as the two lightcone hmits ^ 
and ^ do not commute, the integral 
over k^ of the distribution has a finite de- 
pendence on the regularization parameter 77, 



dk±_ f{x,k±,fi,i]) 



Fix,p,v) (3) 
ordinary pdf , 




which makes the relation with the standard 
operator product expansion not so transpar- 
ent. 

An alternative route is based on the sub- 
tractive method [90l [97]. In this case the 
direction n is kept to be lightlikc but the 
divergences are canceled by multiplicative, 
gauge-invariant countcrterms given by vac- 
uum expectation values of eikonal opera- 
tors. The countcrterms contain in general 
both lightlike and non-lightlike eikonals. For 
this reason studies of u-pdf's in this frame- 
work introduce auxiliary eikonals in direc- 
tion u = (u"*", u~, 0_l) with u"*" and u~ 
nonzero ^71 [57]. The counterterms have 
compact all-order expressions in coordinate 
space. At one loop they provide the reg- 
ularization of the X — > 1 endpoint through 
a kj^ 7^ extension of the plus-distribution 
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regularization, whose specific form is deter- 
mined in [?7] . 

The subtractive method is more system- 
atic than the cut-ofF, and Ukely more suit- 
able for using unintcgratcd parton distribu- 
tions at subleading-log level. In particular, 
the specific choice of the counterterms [27] is 
such that the dependence on the non-lightlikc 
eikonal u cancels in the matrix element at 
y± ~ 0, corresponding to the ordinary pdf. 
The subtractive method may be helpful for 
global NLO analyses incorporating Sudakov 
resummation |98|, I99j . and construction of 
parton-shower algorithms beyond leading or- 
der [U UnOI. See [Ml dni una for stud- 
ies of subtractive methods similar to those 
of [90l [97] in relation to effective-theory tech- 
niques. 

3.3 Further issues in higher order 

Further issues arise in higher order. Non- 
universality of kj_-dcpendcnt parton distribu- 
tions [IH] has recently been studied [^-[5^ 
in the hadroproduction of back-to-back high- 
Pt particles. Potential factorization- breaking 
effects from soft gluons coupling initial and 
final states are shown in Fig. [8] [30] . These 
would enter at a high order of perturbation 
theory (N^LO correction to dihadron pro- 
duction) . Although the role of these correc- 
tions is yet to be fully established, it is inter- 
esting that Coulomb/radiative mixing terms 
also appear to break color coherence [104j in 
high-order contributions to cross sections for 
dijets with a gap in rapidity. A satisfactory 
understanding of factorization should likely 
require a full clarification of these issues. 

4 Conclusion 

Monte-Carlo event generators are being de- 
veloped based on gauge-invariant definition 
of u-pdf's at small x (Sec. 2.1). Compared 
to standard parton shower approaches, these 
Monte-Carlos have the advantage of taking 
into account QCD initial-state radiation pro- 




Figure 8: Soft gluon exchange with spectator 
partons. 

cesses that depend on the largc-k^ tail of par- 
tonic distributions and matrix elements. 

Such processes are likely to be relevant to 
the simulation of complex final states charac- 
terized by multiple hard scales. An example 
is the DIS multi-jet correlations discussed in 
Sec. 2.3. But we expect analogous consid- 
erations to apply for multi-jet final states in 
hadron-hadron collisions at LHC energies. 

k_L Monte-Carlos arc currently being tuned 
and validated using HERA and Tevatron 
data, as discussed in Sees. 2.2 and 2.3. Deter- 
minations of the unintegrated gluon density 
are performed. 

It will be important to extend the notion 
of u-pdf's beyond x <C 1 (Sec. 3) to turn the 
above event generators into general-purpose 
tools. 

This involves a challenging theoretical pro- 
gram, raising new questions on operator ma- 
trix elements (Sec. 3.1), factorization and 
lightconc divergences (Sec. 3.2), and possibly 
new Coulombic effects (Sec. 3.3). 
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